泊松类型方程边界元解法

本文采用高阶拉普拉斯算子基本解将泊松类型方程的区域积分全部变换成边界积分，使计算问题的维数减少一维．通过斯托克斯方程的算例，表明本文所用的方法是有效的方法。     

A Particle-Mesh Method for the Shallow Water Equations Near Geostrophic Balance

In this paper we outline a new particle-mesh method for rapidly rotating shallow water flows based on a set of regularized equations of motion. The time-stepping method uses an operator splitting of the equations into an Eulerian gravity wave part and a Lagrangian advection part. An essential ingredient is the advection of absolute vorticity by means of translated radial basis functions. We show that this implies exact conservation of enstrophy. The method is tested on two model problems based on the qualitative features of the solutions obtained (i.e., dispersion or smoothness of potential vorticity contours) as well as on the increase in mean divergence level.     

The Properties of Voronoi Tessellations for 2D Ising Fluid Near the Second-Order Magnetic Phase Transition

The topological and geometrical properties of Voronoi cells generated for 2D fluid of hard disks with Ising-like spins near the second-order phase transition from the paramagnetic to ferromagnetic phase are described for different disk densities. The comparison with Voronoi cells generated for the random hard disk system is given.     

The effects of flow split ratio and flow rate in manifolds

This paper reports a combined experimental and numerical investigation of three-dimensional steady turbulent flows in inlet manifolds of square cross-section. Predictions and measurements of the flows were carried out using computational fluid dynamics and laser Doppler anemometry techniques respectively. The flow structure was characterized in detail and the effects of flow split ratio and inlet flow rate were studied. These were found to cause significant variations in the size and shape of recirculation regions in the branches, and in the turbulence levels. It was then found that there is a significant difference between the flow rates through different branches. The performance of the code was assessed through a comparison between predictions and measurements. The comparison demonstrates that all important features of the flow are well represented by the predictions.     

Decay of confined,two-dimensional,spatially periodic arrays of vortices: A numerical investigation

The disarrangement of a perturbed lattice of vortices was studied numerically. The basic state is an exponentially decaying, exact solution of the Navier-Stokes equations. Square arrays of vortices with even numbers of vortex cells along each side were perturbed and their evolution was investigated. Whether the energy in the perturbation grows somewhat before it decays or decays monotonically depends on the initial strength of the vortices of the basic state, the extent of lateral confinement and the structure of the perturbation. The critical condition for temporally local instability, i.e. the critical amplitude of the basic state that must be exceeded to allow energy transfer from the basic state to the perturbation, is discussed. In the strongly confined case of a square lattice of four vortices the appearance of enchancement of global rotation is the result of energy transfer from the basic state to a temporally local unstable mode. Energy is transferred from the basic state to larger-scaled structures (inverse cascade) only if the scales of the larger structures are inherently contained in the initial structure of the perturbation. The initial structure of the double array of vortices is not maintained except for a very special form of perturbation. The facts that large scales decay more slowly than small scales and that, when non-linearities are sufficiently strong, energy is transferred from one scale to another explain the differences in the disarrangement process for different initial strengths of the vortices of the basic state. The stronger vortices, i.e. the vortices perturbed in a manner that increases their strength, tend to dominate the weaker vortices. The pairing and subsequent merging (or capture) of vortices of like sense into larger-scale vortices are described in terms of peaks in the evolution of the square root of the palinstrophy divided by the enstrophy.     

潜热型热功能流体强化换热分析

本文建立了分析带有相变微胶囊的潜热型热功能流体的流动和换热过程数理模型,应用有限差分法和移动热源法进行联合求解。计算结果表明,相变微胶囊的加入:较好地提高了流体的换热性能。获得了相变颗粒大小和体积分数对强化换热影响等结果。计算结果对该类流体的设计和应用提供了理论依据。     

The modified characteristic finite difference fractional steps method for the coupled system of fluid dynamics in porous media and its analysis

For the coupled system of multilayer fluid dynamics in porous media, the modified characteristic finite difference fractional steps method applicable to parallel arithmetic is put forward and two‐dimensional and three‐dimensional schemes are used to form a complete set. Some techniques, such as calculus of variations, energy method, piecewise biquadratic interpolation, multiplicative commutation rule of difference operators, decomposition of high order difference operators and prior estimates are adopted. Optimal order estimates in L² norm are derived to determine the error in the approximate solution. This method has already been applied to the numerical simulation of multilayer fluid dynamics in porous media. © 2003 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 19: 665–681, 2003.     

A new VOF‐based numerical scheme for the simulation of fluid flow with free surface. Part II: application to the cavity filling and sloshing problems

Finite element analysis of fluid flow with moving free surface has been performed in 2‐D and 3‐D. The new VOF‐based numerical algorithm that has been proposed by the present authors (Int. J. Numer. Meth. Fluids, submitted) was applied to several 2‐D and 3‐D free surface flow problems. The proposed free surface tracking scheme is based on two numerical tools; the orientation vector to represent the free surface orientation in each cell and the baby‐cell to determine the fluid volume flux at each cell boundary. The proposed numerical algorithm has been applied to 2‐D and 3‐D cavity filling and sloshing problems in order to demonstrate the versatility and effectiveness of the scheme. The proposed numerical algorithm resolved successfully the free surfaces interacting with each other. The simulated results demonstrated applicability of the proposed numerical algorithm to the practical problems of large free surface motion. It has been also demonstrated that the proposed free surface tracking scheme can be easily implemented in any irregular non‐uniform grid systems and can be extended to 3‐D free surface flow problems without additional efforts. Copyright © 2003 John Wiley & Sons, Ltd.     

The dynamics of a two-component fluid in the presence of capillary forces

In the present paper we study the qualitative behavior ast→∞ of the solution of the Cauchy problem for a system of equations describing a dynamics of a two-component viscous fluid. The model under consideration takes into account the mutual diffusion of the fluid components as well as their capillary interaction. We describe the ω-limit set of trajectories of the dynamical system generated by the problem. It is proved that the stationary solution of the problem, is a homogeneous stationary distribution of one of the components, is asymptotically stable. Any other stationary solution is not asymptotically stable and is even unstable if there are no close stationary solutions corresponding to a smaller energy level. Translated fromMatematicheskie Zametki, Vol. 62, No. 2, pp. 293–305, August, 1997. Translated by A. M. Chebotarev